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Abstract 

In 1952 Bohm presented a theory about non-relativistic point-particles that move 
deterministically along trajectories and showed how it reproduces the predic- 
tions of standard quantum theory. This theory was actually presented before 
by de Broglie in 1926, but Bohm's particular formulation of the theory inspired 
Epstein to come up with a different trajectory model. The aim of this paper is to 
examine the empirical predictions of this model. It is found that the trajectories in 
this model are in general very different from those in the de Broglie-Bohm theory. 
In certain cases they even seem bizarre and rather unphysical. Nevertheless, it is 
argued that the model seems to reproduce the predictions of standard quantum 
theory (just as the de Broglie-Bohm theory). 

1 Introduction 

In the de Broglie-Bohm theory, discovered by de Broglie [1] and rediscovered (in a differ- 
ent form) by Bohm [2,3], an individual closed system is described by its wavefunction 
ijj(x, t) on configuration space, with x = (xi, . . . , Xjv), which satisfies the non-relativistic 
Schrodinger equation 

ifc^iM = _f- ^-V&Ks, t) + V(x)^(x, t) , (1) 

k=l 

and by particle positions xi(t), . . . , Xjv(t), whose possible trajectories are solutions of 
the guidance equations 

dxi. 1 „ 

-± = —V k S, 2 
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where S is the phase of the wavef unction, that is, = exp(iS/H). 

An important property of the guidance equations is that they preserve the distri- 
bution \ip\ 2 . That is, if the distribution of particles is given by \ip(x,t )\ 2 at a certain 
time to, then the distribution is given by \ip(x, i)\ 2 at other times t. This property is 
called equivariance [4] (and \ip\ 2 is actually the only distribution that is a suitably local 
functional of ip that satisfies it [5]). Equivariance essentially follows from the fact that 
the distribution \if)(x, t)\ 2 satisfies the continuity equation 



just as any distribution p(x,t) that is transported along the de Broglie-Bohm trajecto- 
ries. The distribution \ip\ 2 plays the role of an equilibrium distribution (see for example 
[4,6,7]) and is called the quantum equilibrium distribution. Given the quantum equi- 
librium distribution and the fact that measurement results are generally ultimately 
recorded in positions of macroscopic pointers, like instrument needles, computer out- 
print, etc., it almost follows immediately that the de Broglie-Bohm theory reproduces 
the standard quantum mechanical predictions (see for example [3,8-11]). 

Instead of following the above presentation of the theory, which is along the lines of 
that of de Broglie, Bohm actually presented the theory in a Newtonian form. Bohm ob- 
served that by using the polar decomposition of the wavefunction, i.e. ifi = exp(iS/h), 
the Schrodinger equation yields an equation that is similar to the Hamilton- Jacobi equa- 
tion of Newtonian mechanics. This inspired Bohm to adopt a Newtonian type equation 
with an extra -^-dependent potential as the basic equation of motion, thereby regarding 
the guidance equations (2) as constraints on the possible momenta. 

It was this alternative formulation of Bohm that inspired Epstein to come up with 
an alternative trajectory model [12, 13]. Epstein found this model by considering the 
wavefunction in momentum space rather than in configuration space and by following 
an analysis similar to that of Bohm. 

The aim of this paper is to examine the empirical predictions of Epstein's theory. 
To our knowledge this has not been done before. It is found that localized macroscopic 
systems, such as pointers in measurement situations, seem to be located at locations 
predicted by standard quantum theory and this with the same probabilities as given 
by standard quantum theory. This indicates that Epstein's trajectory theory seems to 
reproduce the standard quantum predictions, just as the theory of de Broglie and Bohm. 

2 Epstein's trajectory model 

Just as Bohm, Epstein derived his model by considering the formal analogy with classical 
Hamilton- Jacobi theory. In presenting Epstein's model we will not consider this analogy, 
but instead restrict ourselves to give just the basic equations of motion (just as (1) and 
(2) can be regarded as the basic equations of motion in the de Broglie-Bohm theory). 




(3) 
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2.1 Equations of motion 



In Epstein's trajectory model the wavefunction is written in the momentum representa- 
tion, by taking the Fourier transform of the wavefunction in the position representation: 



if>(p,t) = F p {ip(x,t)) = (2nh)-^ / dxe- ix - p/h ^(x,t) . 



(4) 



Application of the Fourier transform to both sides of the non-relativistic Schrodinger 
equation (1) yields 



N 



2 



k=i 



2m k 



This implies 



where 



dt 



+ I+(p,t) = 0, 



rffat) = ^Re (ip(p,t)T p (V(x)il>(x,t)j) . 



(5) 



(6) 



(7) 



For Epstein's model the quantity 1^ needs to be written as a divergence J2k ^fc'Jfc > where 
V k = {d/dp k x,d/dp ky ,d/dpkz) and p fc = {pkx,Pk y ,Pkz), in such a way that the above 
equation (6) can be interpreted as a continuity equation for the momentum density j^) 2 
with currents jjj? (in particular the currents should die off sufficiently fast for \p\ — > oo so 
that there is no flux at infinity). The possibility of introducing such currents depends 
of course on the form of the potential. This is discussed below. 

Given currents j]?, the first step in Epstein's model is to introduce the variable 3 
p = (pi, . . . , ptv) whose time evolution is determined by the differential equations 



dpk 

dt 



(8) 



The next step is to introduce particle positions xi, . . . , x^, whose trajectories are defined 
as 

'r(p,t)(ihV k )^(p,t) 



x fc (i) = Re 



-V k S(p,t) 



p=p(t) 



p=p(t) 



(9) 



where the polar decomposition ip = \ exp(iS/h) is used and where the expressions on 
the right hand side are evaluated for an actual trajectory of the variable p. 4 



3 One could call these variables momenta. However, it should be kept in mind that the relation 
Pfc = rnxfc, with Xfe the fc-th particle position as defined in (9), in general does not hold. The exact 
relation to the velocity could be found by differentiating the expression in (9) with respect to time. 

4 These expressions correspond to Holland's local expectation values [9] for the position operators 
Xfc in the momentum representation, where x^ — > ifiV^ , evaluated for the actual value of p. 
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So in Epstein's trajectory model systems are described by the triplet (x,p, ip) (where 
the particle positions constitute the primitive ontology [4,14]). The wavefunction tp 
determines the time evolution of the variable p and the wavefunction together with the 
variable p determine the configuration x. 

Let us now return to the construction of the currents j^. The construction of such 
currents was discussed in detail in [15]. If the potential is given by a differential operator 
in momentum space, that is, if there exist functions V n (p) such that 

F p (V(x)i;(x, t)) = J2 V n (p)D n ^(p, t) , (10) 

n>0 

where n = (n lx , n ly , n lz , n Nx , n Ny , n Nz ) e N m and D n = (^/dp^, d n ^/dp^), 
then currents can be found of the form 

3tM= E h,nm(p)D n ^( P ,t)D m r(p,t), (11) 
n,m>0 

where Jk,nm are certain functions that depend on the V n . 

In the case the potential is not given by a differential operator in momentum space, 
for example if it contains a Coulomb potential, then one could solve the Poisson equation 
V 2 F = 1^, where V 2 = Xlfc^fc * s the 3iV-dimensional Laplacian, and define = V^F. 
A solution is given by F = V -2 /^ [15] (given that it is well-defined), in which case 

tf(P,*) = V*ii*(p,f). (12) 

This method of constructing a current was suggested by Epstein himself [13] in response 
to a worry raised by Bohm [16] concerning the treatment of certain potentials such as 
the Coulomb potential. In the following we just assume that there exists a suitable 
choice of current, so that the dynamics (8) for the variable p can be introduced. 



2.2 Ensemble distribution 

It is assumed that the distribution of the variable p over an ensemble of systems all 
described by the same wavefunction ip is given by \ip\ 2 . This distribution is equivariant 
with respect to the time evolution of the variable p. 

The implied distribution for the configuration x is then given by 

P+{x,t) = J dpm P ,t)\ 2 5[x + VS(p,t)}, (13) 

where we defined 

N 

6[x + VS(p, t)} = J] 5[x fe + V k S(p, t)} . (14) 
k=i 

The distribution could be regarded as an equilibrium distribution, analogous to the 
quantum equilibrium distribution |-?/>| 2 in the theory of deBroglie and Bohm. However, 
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the distribution p^ is in general very different from \ip\ 2 (this is illustrated in an example 
below). Only the latter equals the quantum mechanical position distribution. 

While the position distribution p^ is in general different from the quantum mechan- 
ical one, the expectation value for position equals the quantum mechanical one. 5 That 
is, if we write x fc = (x k ,y k ,z k ), then 



(xk) = j dxp^(x)x k = - j dp\i)(p)\ 2 -^-(p) = J dpip*(p) ^ih 



dp 



kx 



(15) 

and similarly for the other Cartesian coordinates. We also have that the standard 
deviation satisfies the inequality Ax k ^ Axfc, and similarly for the other Cartesian 
coordinates. This follows from 

(4) = (4) -h 2 jd P (^H^)l) 2 < (4) . (16) 

2.3 Example 

Consider a free particle as an example. Its wavefunction in momentum space is given 

iKp,*) = f\v)e- ipt/2mh (17) 

where / is an arbitrary (square integrable) function. Suppose further that / is real. 
Since the potential is zero the evolution equation (8) reduces to dp/dt = 0, so that the 
possible trajectories are given by 

*(«) = * (18) 

where p is a constant vector, distributed according to |^(p,t)| 2 = f{p) 2 over an en- 
semble. Note that all the trajectories pass through the origin at t — 0. So clearly the 
initial position distribution p^(x, 0) = 5{x) is different from the quantum mechanical 
position distribution |^(x, 0)| 2 . Also the trajectories are significantly different from the 
trajectories in the theory of de Broglie and Bohm. For example the trajectories in the 
latter model never cross in configuration space, while our example illustrates that the 
trajectories in Epstein's model are clearly allowed to do that. 

The particular position distribution in Epstein's model seems to lead to rather bizarre 
consequences. Consider again a free particle, with wavefunction in the superposition 

0(x, t) = ^(^(x - a, t) + ^(x + a, t)) , (19) 

where ^(x, t) is a wavefunction localized near the origin at time t — and where a is 
some constant vector. In the momentum representation this becomes 

0(p, t) = V2 cos(a • p//#(p, t) . (20) 



'This is a defining property of Holland's local expectation value. 
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With ip of the form (17) with / real, the particle distribution for an ensemble is p^(x, 0) = 
5(x). While we would have expected that the particles are typically located around a 
or —a at time t — 0, they are located at the origin, independently of the value of a. For 
example, the wavefunction of the particle in the double slit experiment is of the form 
(19) (after it has passed the slits) and the above analysis implies that the actual particle 
position would initially come from the midpoint between the slits. Hence, Epstein's 
model seems to yield a bizarre and rather unphysical picture for isolated particles in 
certain cases. Nevertheless we will argue in the next section that the model actually 
seems capable of reproducing the predictions of standard quantum theory, because it 
can account for positions of macroscopic systems. 



3 Reproducing the predictions of standard quantum 
theory 

3.1 Macroscopic systems 

Consider a wavefunction ip(x) that represents a macroscopic system, like for example an 
instrument needle or a cat. This wavefunction has most of its support on configurations 
that on the macroscopic level give the image of the macroscopic system. That is, those 
configurations all correspond to the same macrostate of the system. Let us call this set 
of configurations B. In the theory of deBroglie and Bohm, where the configurations are 
distributed according to (^(a;)! 2 over an ensemble, the configurations are typically in B. 
So the particles typically yield the image of the macroscopic system. This is true also 
for Epstein's theory: configurations are typically within B because the distribution 
satisfies (xk) = (xk) and Axk ^ Ax fc , cf. Section 2.2, and because B can be assumed to 
have a suitably simple form. 

Consider now wavefunctions that are in a superposition of macroscopically different 
states of the macroscopic system. First assume that the wavefunction ip(x) represents 
a macroscopic system that is localized near the origin of our reference frame in phys- 
ical 3-space and call B again the set of configurations that on the macroscopic level 
yield the image of the macroscopic system. The wavefunction ip a (x) = ip(x — a), with 
a = (a, . . . , a), will represent the same system located near a. Consider now the super- 
position 

N 

<P(x) = ^=(Mx)+^-a(x)), (21) 

where N is a normalization factor. If |a| is sufficiently large, the wavefunctions ip a and 
i[)- a will only have negligible overlap, so that N ~ 1. This is definitely the case if |a| 
is much bigger then the size of the macroscopic system. In the theory of de Broglie 
and Bohm the fact that those wavefunctions have only negligible overlap implies that 
the configuration will typically display the macroscopic system localized at either a or 
—a. This is how the theory deals with Schrodinger's cat paradox. If the wavefunction 
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of a system is in a superposition of macroscopically distinct states, then the particle 
configuration will typically display only one of those macroscopic states. 

In Epstein's theory the situation is a bit different. Let us analyze what the typical 
particle configuration is for the state 0. First of all, the Fourier transform of is given 

(j){p) = y/2N cos(a • p/h)ip(p) , (22) 
where ip is the Fourier transform of ip. Hence the position distribution is given by 

p*(x) = 2N 2 J dp cos 2 (a • p/K) \i){p) \ 2 5[x + VS(p)] , (23) 

where S is the phase of ip. Since 

p*(x) ^ 2N 2 J dp\i>(p)\ 2 5[x + VS(p)] = 2N 2 p^(x) (24) 

and since p^ has negligible support outside B (and N 1), also p^ must have negligible 
support outside B. Hence the position configuration typically does not display the 
macroscopic systems localized near a or —a, but rather near the origin. 

This feature seems to suggest that Epstein's model will fail to reproduce the pre- 
dictions of standard quantum theory. However this is not so. The reason is that we 
did not include the description of the environment in our analysis. The introduction 
of the environment drastically changes the description. To see this, consider a general 
superposition 

4>'(x,x e ) = Cit/Ji(x)xi(x e ) + c 2 ^ 2 (x)x 2 (x e ) , (25) 

where the Cj, i — 1, 2, are constants, where the normalized states ipi are macroscopically 
different states of a localized macroscopic system (for example ipi = ip a and ip2 = V'-o) 
and where the normalized states Xi represent that part of the environment that is in 
direct interaction with the macroscopic system, like for example the air molecules that 
scatter off the system. This implies that the Fourier transforms Xi an d X.2 have negligible 
overlap in momentum space. It is namely sufficient that there is one air molecule in the 
environment that scatters off the macroscopic system and obtains a different momentum 
as a result. In a realistic situation there will be many molecules that obtain a different 
momentum. As a result, the momentum distribution is approximately given by 

l0We)| 2 « |c i r|^( P )| 2 |Xl(Pe)| 2 + \c2\ 2 \M P )\ 2 \MPe)\ 2 , (26) 

with |ci| 2 + |c 2 | 2 ~ 1 (when the Xi are completely non-overlapping these relations hold 
exactly). The corresponding position distribution is then given by 

~ J dpdp e (^ Cl \ 2 \^M\ 2 S\x + VSM}\Xi^ (27) 

+ \C2\ 2 \MP)\ 2 S[X + VS 2 (p)]\x2(pe)\ 2 6[x e + V e S e ,2( Pe )}^ (28) 

= | Cl | 2 p^(x)p^(x e ) + \c 2 \ 2 p^(x)p^(x e ) , (29) 
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where S e j are the phases of Xii i = 1? 2 and V e represents the derivative with respect 
to p e . Since the densities p^ 1 and p^ 2 have approximately different supports (they have 
there support approximately concentrated on different macrostates), the configuration 
of the macroscopic system will typically be in the macrostate corresponding to either 
■01 or ip2 and this with probabilities respectively given by |ci| 2 and | C2 1 2 , which are 
the same probabilities as in standard quantum theory. For example, turning back to 
the particular superposition (21), the above analysis shows that when the environment 
is taken into account, the configurations will typically display the macroscopic system 
localized either near a or —a, with equal probability. 

Note that in many cases, the states ip\ and ip 2 are themselves non-overlapping. It 
is namely sufficient that one particle of the system has a different momentum in ipi 
compared to ip 2 (this will for example definitely be the case for a superposition of a live 
and dead cat). In such a case there is actually no need to invoke the environment to 
establish that the particle configuration will typically display one of the macroscopic 
states. 

3.2 Collapse of the wavefunction 

In the model of Epstein, the wavefunction evolves according to Schrodinger's equation 
at all times. It never collapses. However, under certain circumstances the wavefunction 
can undergo an effective collapse. Effective collapse means the following. Suppose the 
wavefunction is given by a superposition ip = ipi + ip 2 . Then it can happen that the 
dynamics of the actual configuration depends only on the wavefunction ipi or ip 2 from 
some time onwards. Since the other wavefunction then does not play any further role 
in the dynamics of the configuration, it can be removed from the description of the 
evolution of the configuration. This is called an effective collapse. 6 

The circumstances under which there is effective collapse depend on the form of 
the dynamics. In Section 2, two possible types of dynamics were considered for the 
variable p. For a potential that is given by a differential operator in momentum space, 
the current (11) could be chosen for the equation of motion of the variable p. For a 
potential that is not of such a form the current (12) could be chosen (given that the 
expression is well-defined). 

In the case the potential is a differential operator, there is effective collapse if the 
wavefunction is given by a superposition 

$(p,t)=$i(p,t)+$ 2 (p,t), (30) 

where ipi and ip 2 have negligible overlap in momentum space from some time t onwards. 

6 In the de Broglic-Bohm theory, the wavefunction of a subsystem of the universe can be defined as 
the conditional wavefunction [4] (which is constructed from the universal wavefunction and the actual 
positions of the particles not belonging to the subsystem). This wavefunction actually undergoes what 
could be called an actual collapse [4] . In the context of Epstein's theory, the wavefunction of a subsystem 
could be introduced in a similar way (using the universal wavefunction in momentum space and the 
actual p values not belonging to the subsystem). However, we will not pursue this further here. 



8 



This is because the current (11) then has a similar decomposition: 

f(p,t)~j^(p,t)+j^(p,t), (31) 

from time to onwards. This implies that if the momentum variable p is within the 
support of ipi, I = 1,2, at time to, then its future j^volution will be determined only by 
ipi, through the equation of motion dp/dt = j^ l /\ipi\ 2 - In addition we have that 

K k (t)^-^kSi{p,t)\ p=p(t) (32) 

with Si the phase of ipi. Hence the evolution of the particles only depends on the 
wavefunction ipi, so that the other packet in the superposition can be removed from the 
description of the evolution of the particles. Since the variable p is distributed according 
to \ip\ 2 over ensembles, it is easily verified that the probability with which an effective 
collapse to ipi occurs is the same as the probability standard quantum theory would 
assign for an actual collapse of the wavefunction to tpi. 

In the case the potential is not a differential operator in momentum space and the 
current (12) is used in the equations of motion for the variable p, it is not so clear 
whether the above conditions on the state (30) are sufficient to guarantee effective 
collapse. If, in addition to having negligible overlap, the two wavefunctions ipi both 
satisfy the Schrodinger equation (in the case the potential is given by a differential 
operator the latter automatically follows from the non- overlap), then I* I^ 1 + J^ 2 . 
From the equations d\ipi\ 2 /dt + I^ 1 = it follows also that the I^ 1 have approximately 
the same support as ipi (provided ipi does not fluctuate too fast outside the region 
where most of its support is concentrated). As such the functions I^ 1 and J^ 2 will have 
negligible overlap. However the currents (12) are non-local in momentum space so that 
the latter feature does not immediately imply that the components j^ 1 and j^ 2 in the 
current ~ j^ 1 + j^ 2 are non-overlapping. While the currents j^ 1 die off outside the 
region where most of the support of I^ 1 is concentrated (the currents j^ 1 relate to I^ 1 in 
the same way the electric field relates to the charge density by which it is generated), a 
more detailed analysis would be required to find the rate at which this happens. While 
equation (32) holds regardless of this, the configuration will depend on both terms of the 
wavefunction through the actual value of p, if the time evolution of the latter depends 
on both terms. 

3.3 Measurement situation 

Consider now the description of an ideal quantum measurement in terms of Epstein's 
theory. Before the measurement it can be assumed that the wavefunction tp is given by 
a product wavefunction 

i> = i>si>mi>e , (33) 

where the wavefunction ip s represents the microscopic system under observation, where 
ip m represents the measurement device which includes some macroscopic pointer, for 
example an instrument needle, and where ip e represents that part of the environment 
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that is in direct interaction with the measurement device. During the measurement the 
wavefunction will evolve into a superposition, i.e. 



where the wavefunctions ip 8j i represent the different eigenstates of the quantum oper- 
ator that is being measured, where the ip m j represent the corresponding states of the 
measurement device, including the macroscopic pointer indicating the measurement out- 
come, where the ip e> i represent the different associated states of the environment, and 
where the q are constants. 

As explained in Section 3.1, it is guaranteed that the different packets ipi = ipsjipmjipej 
have negligible overlap in momentum space, since at least the packets i/j e j have negligible 
overlap. And as a result, the configuration x m of the measurement device will typically 
display one of the possible states of the macroscopic pointer, thereby displaying the 
outcome of the measurement. The possible results are the same as those expected ac- 
cording to standard quantum theory and are obtained with the same probabilities as 
predicted by standard quantum theory. As such the standard quantum predictions for 
measurements are reproduced by Epstein's theory. 

In the case the potential is given by a differential operator in momentum space, the 
non-overlap of the packets ipi guarantees that there will be an effective collapse to one of 
these packets with the same probability as in standard quantum theory. In this way the 
collapse postulate of standard quantum theory is recovered. In the case the potential is 
not given by a differential operator, it is unclear whether there is an effective collapse. 
If there is no effective collapse it might be that the actual configuration moves (very 
fast) between the configurations that display the various macroscopic states. While 
this feature would be empirically unverifiable because the configurations of the system, 
the measurement device and the observer's memory would maintain their correlation, 
it seems unacceptable for a physical theory. In any case, even if there is no effective 
collapse, the model reproduces the predictions of standard quantum theory because, 
again, the possible measurement results are the same as those in standard quantum 
theory and they are obtained with the same probabilities. 

4 Conclusion 

While Epstein's trajectory model is very different from that of de Broglie and Bohm, it 
seems to agree with the latter on the empirical level. Examples of other such trajectory 
models have been given before [17-19]. Epstein's model differs at least in one aspect 
from those theories, namely that the position distribution is in general very different 
from \ip\ 2 . An example of another theory which shares this property is that of de 
Polavieja [20,21], which is developed using the Schrodinger equation in a phase space 
representation. The extent of agreement of the latter theory with standard quantum 
theory still needs to be investigated. 




(34) 
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Similar models could be constructed starting from representations intermediate to 
the position and momentum representation. After introducing a variable in such an in- 
termediate representation, the position variable would probably needed to be introduced 
through some construction like Holland's local expectation value. Of all these theories, 
the one of deBroglie and Bohm is arguably the simplest and most natural. It doesn't 
require a variable intermediate between the wavefunction and the position variable (like 
the momentum variable in Epstein's theory). 
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